Raman Model Predicting Hardness of Covalent Crystals 
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Abstract 

Based on the fact that both hardness and vibrational Raman spectrum depend on the intrinsic property 
of chemical bonds, we propose a new theoretical model for predicting hardness of a covalent crystal. The 
quantitative relationship between hardness and vibrational Raman frequencies deduced from the typical 
zincblende covalent crystals is validated to be also applicable for the complex multicomponent crystals. This 
model enables us to nondestructively and indirectly characterize the hardness of novel superhard materials 
synthesized under ultra- high pressure condition with the in situ Raman spectrum measurement. 
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Design and synthesis of new superhard materials are of great interest due to their numerous 
appKcatioa, fl^flfl Hardness, as an „n port an t _pic pn^ica! pr ope rty , „ nn^ood 
as the resistance offered by a given material to applied mechanical action. For a crystalline material, 
hardness is an intrinsic property. Its prediction, from the microscopic electronic structure, is a 
crucial issue and a powerful challenge in condensed matter physics and materials science. Recently, 
the microscopic model connecting hardness with the nature of chemical bond has shed light on 
the quantitative estimations of other macroscopic properties. |(], |?], To predict hardness using 
the above models, the exact crystal structure must be known. Experimentally, the ultra-high- 
pressure tech niq ue such as diamond anvil cell (DAC) is extensively used to explore new superhard 



materials. 



3. [9, 10, [ill . [12I Due to the limited size of the sample synthesized in DAC, it is often 
very difficult to either conclusively determine the exact atomic arrangement or the hardness for 
some newly synthesized samples. These quantitative models are thereby limited on some practical 
applications. However, the in situ measurements of Raman and infrared equipped on DAC can 
provide the diagnostic electronic features of the phase. 

Raman spectroscopy is commonly used in chemistry, since vibrational information is specific 
for the chemical bonds in a crystal. It therefore provides a fingerprint for identifying the crystal. 
The Raman scattering process involves the change in susceptibility while the infrared absorption 
connects with the change in dipole moment. The Raman modes include vibrational, rotational, 
and librational modes, while the infrared absorption mode involves translational mode only. Fur- 
thermore, the vibrational Raman modes can also be classified into the longitudinally optical (LO) 
and transversely optical (TO) modes, originating from the tensing/compressing and bending of 
chemical bonds, respectively. [5[j Based on the microscopic understanding, hardness is naturally the 
resistance of chemical bond per unit area to indenter.[6] In the experimental measurement by using 
the indenter, the hardness value is determined from the ratio of the load to the indentation area. [5] 
This resistance is related to the tensing/compressing and bending, but not the rotating, librating 
or translating of chemical bonds because the bond length will keep the same for the translation and 
rotation (libration) modes, it wouldn't embody the effective bond strength if considering simulation 

of the indentation, the contacted bonds would be deformed or broken. 

' 1 — I 

In the previous works, the resistant force of bond can be characterized by energy gap, [(J ref- 
erence energy [?[] and bond electronegativity. [8] In our understanding, the vibrational Raman fre- 
quency/energy can also be used to describe the resistant force, due to the reasons as mentioned 
above. In this work, we create a quantitative relationship between hardness and vibrational Raman 
modes in nonresonant first-order Stokes Raman spectrum. Raman spectra find other important 
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application on the prediction of hardness. 

Inasmuch as the vibrational Raman frequency of a chemical bond increases as the force constant 
of a chemical bond increases, the vibrational Raman frequency embodies is in fact the bond strength 
of a chemical qualitatively. The vibrational Raman frequency can be determined by the first-order 



aecause it lies on the frequency u of the optical phonon at the Brillouin-zone center 
Raman intensity of the Raman mode with the Raman frequency u) can be 



response, 
(r point) 

calculated under the Placzek approximation. [141] Raman spectra can be calculated by the PWSCF 
implementation within the density functional perturbation theory (DFPT) formalism. [141] 

In the previous works, twice of the band gap or band energy, [g] the reference energy, [3] and 
electron-holding energy[8[ were used to define the resistant force of bond. We introduce the Raman 
frequency of a vibrational Raman active bond, u m , to embody the resistant force F m of this bond 
as 

F m oc w m exp(w m /u;o), (1) 

where ojq is a constant. The suffix m indicates the ordinal number of the vibrational Raman 
mode. Inasmuch as hardness of a material is attributed to the collective contribution of all the 
chemical bonds in any direction, hardness can be described by the compromise resistant force. 
Due to differences in the resistant force among different types of bonds, the compromise resistant 
force provided by all the chemical bonds with the vibrational Raman active property should be 
depicted by a weighted geometric average [6] of all the vibrational Raman frequencies. Thus the 
contributions of all the chemical bonds can be equivalent to that of a single isotropic chemical bond 
with an imaginary Raman frequency Ud, here we call oJd the diagnostic Raman frequency, which is 
expressed as 



Here I m is the relative intensity of the mth vibrational Raman mode with the Raman frequency 
uj m and indicates its contribution fraction or the weighted factor. Therefore, the single equivalent 
isotropic chemical bond provides the resistant force F should be 

F oc Lddex.p(uj d /uJo). (3) 

It is clear that F — > when iv<i — * 0. Because hardness is understood as the resistance offered by a 
given material to applied mechanical action, the Vickers hardness can be expressed as Hy = AF, 
where A is a proportional constant. Despite the fact that we cannot give the ab initio deduction 



for the exact expressions of A and ujq from the theory, we can still obtain their values based on 
the semi-empirical method. It is well known that the zincblende crystals have the simplest Raman 
spectra with two Raman modes (one LO and one TO) and the zincblende diamond is the hardest 
material so far. Therefore, we can use the known Raman spectra and hardness for the typical 
zinc-blende covalent crystals, to determine the constants of A and ujq. From Hy against for the 
first 14 typical zincblende covalent crystals from Table I as shown in Fig. 1, we obtain 

ify(GPa) = 0.011a; (i exp(u; (i /704.8), (4) 

where ojd is in the unit of cm" 1 . For all the 22 crystals listed in Table I, the Vickers hardness 
values calculated by the above formula are given. Our results are in good agreement with the 
experimental data and the theoretical values by the other models. 

To confirm the universality of Eq. (4) , we first apply it to the wurtzite crystals with two types of 
chemical bonds. The results validate the good agreements between the calculated and experimental 
values as shown by circles in Fig. 2. 

Next we focus on some typical complex crystals as listed in Table II. For /3-SisN4 (a fourfold 

coordinated Si atom is linked by four threefold coordinated N atoms), from the seven vibrational 

□ 

Raman modes measured in experiment. [191] we obtain to^ = 818.9 cm 1 by Eq. (2) and then 
Hy = 28.8 GPa by Eq. (4). The calculated Vickers hardness is in good agreement with the 
experimental value of 30 GPa and the theoretical value of 30.3 GPa.fg] For the experimental 
Raman spectra of a-Si02 in Ref. 14 and Stishovite in Ref. 20, the theoretical analyses reveal that 
the strong Raman mode at 207 cm" 1 for Q-Si02 and the two strong Raman modes at 589 and 
231 cm" 1 for Stishovite should be excluded for evaluating hardness, because they belong to the 
rotational Raman modes. When other very weak Raman modes are ignored simultaneously, only 
two relatively strong Raman modes need to be included for calculating hardness: 464 and 450 cm" 1 
for a-Si02 in Ref. 14 as well as 967 and 753 cm -1 for Stishovite in Ref. 20. The calculated values 
are 9.8 GPa for a-Si02 and 25.3 GPa for Stishovite, respectively, which are in good agreement 
with the experimental values of 11 GPa in Ref. 6 for a-Si02 and of 17-23 GPa in Ref. 5 (but 
lower than 33 GPa in Refs. 5 and 6) for Stishovite. For a-Si02, the Vickers hardness of 30.6 GPa 
predicted by the microscopic model is remarkably higher than the experimental value of 11 GPa 
in Ref. 6, whereas our predicated value of 9.8 GPa is very close to the experimental value. 

We now turn to explore the boron-rich systems, includin g q -Boron, B4C and P^O. Based on 



both the theoretical and experimental analyses, 
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22, 



23 



25] Vast et al. have pointed out that 



the Raman modes of a-Boron could be classified into three groups: intericosahedral modes at high 
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frequency, intraicosahedral modes at middle frequency, and librational modes at low frequency. [22 ] 
This opinion is also feasible to B4C and BgO. Lazzari et al. have validated the two measured 
low-frequency Raman modes (498 and 534 cm -1 ) for B4C are the rotational and librational modes, 



respectively. 



24i | For a-Boron, from the experimental Raman spectrum, [21[ we evaluate the 



Vickers hardness to be 39.8 GPa (oj<i = 945.8 cm 1 ), which agrees with the experimental value of 



42 GPa. [25J By using our present model combined with the measured Raman spectra, 2l|, |23|] the 
calculated hardness are 40.2 GPa (u;^ = 950.9 cm" l ) for B 4 C and 44.9 GPa (uo d = 994.6 cm" 1 ) for 
BgO, which are in good agreement with the experimental values of 42-49 GPa for B4C in Ref. 5 
and 45 GPa for BgO in Ref. 25, respectively. 

Finally, we concern the ter nary superhard BC2N. Recently, we reported a new phase (Z-BC2N) 
with the P-42M space group, [36|] and the simulated XRD pattern is in good agreement with the 
experimental data. [ill] However, the first-principles calculation reveals that the non-resonant first- 
order Stokes Raman spectrum of Z-BC2N does not match to the measured Raman pattern of the 
synthesized BC2N sample. [3] We construct a modified structure of BC2N, named as z*-BC2N 
here, with the same space group as Z-BC2N. In fact, the unique difference between z*-BC2N and 
Z-BC2N is that only B and N atoms are interchanged each other. The simulated Raman spectrum 
of z*-BC2N is in good agreement with the experimental Raman pattern of the synthesized BC2N 
sample. [lj| Based on the simulated Raman spectrum of z*-BC2N as Table II, the evaluated 



m 



hardness of z*-BC2N is 75.5 GPa, which agrees with the experimental value of 76 GPa. 

It should be noted that the role of d valence electron in chemical bond are not considered 
here. This model could also be simplified for roughly estimating hardness of potential superhard 
materials with the Vickers hardness above 40 GPa, as follows: to avoid the procedure eliminating 
the translational and rotational Raman modes, we only need to select the strong Raman modes 
with their Raman frequencies higher than the critical frequency u> c = 687 cm -1 (u c is the diagnostic 
Raman frequency of Hy = 20 GPa, when a substance with its Vickers hardness between 20 and 
40 GPa is usually considered as a hard material). For example, for superhard diamond-like BC5 



(d-BC 5 ).. 



371 ] there is one strongest Raman mode located at 1200 cm above uj c in the measured 



Raman spectrum. By using this single strongest Raman mode, we estimate its Vickers hardness to 
be 72.4 GPa, which is in excellent agreement with the experimental value of 71 GPa. (37}] To display 
an intuitive comparison, Fig. 2 summaries the calculated and experimental Vickers hardness values 
for 25 crystals. The good agreement validates the predictive power of our model. 

In conclusion, we propose a model that reveals a quantitative relationship between hardness 
and vibrational Raman frequencies. This model has a tremendous advantage that, to evaluate the 
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hardness synthesized materials, we do not need to know the exact atom arrangement provided 
that the Raman spectrum can be measured. The present work validates a universal technique 
for the nondestructive and indirect hardness measurement, it also be potential used to explain 
enhancement of the surface or vacancy hardness because of the high sensitivity of the Raman 
spectra to the small change in structure. Moreover, since different Raman configurations for a 
single crystal will result in different Raman spectra, our model can be anticipated to explore the 
anisotropy of the hardness to some extent. This Raman model finds in fact a new application for 
the Raman spectroscopy. 

This work is in part supported by the National Natural Science Foundation of China un- 
der Grants Nos. 50821001 and 50532020, by the 973 Program of China under Grant Nos. 
2006CB921805 and 2005CB724400. 
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Table I. Hardness, LO m and I m (in parenthesis) of LO and TO Raman modes, for 22 kinds of 
diamond-like crystals. Hy ur , Hy ao and Hy m are the calculated Vickers hardness by our, Gao's [6] 
and Simunek's [7] models, respectively. Hy Xp is the experimental Vickers hardness (unless noted, 
from Refs. [6,7]). Despite the fact that there are no measured Vickers hardness for 8 materials 
from ZnSe and below, their Knoop hardness values marked by the asterisk are listed as a reference. 



Crystals 




WTO 


H Our 


H Gao 


H s im 


rjExp 

n v 


Diamond 1332 (1.0) a 1332 (1.0) a 


97.0 


93.6 


95.4 


96 


BN 


1305 (0.7) b 


1055 (1.0) b 


64.9 


64.5 


63.2 


66 


SiC 


976 (0.6) c 


796 (1.0) c 


32.0 


30.3 


31.1 


34 


BP 


829 (1.0) b 


799 (0.2) b 


29.2 


31.2 


26.0 


33 


A1N 


893 (0.2) b 


668 (1.0) b 


20.9 


21.7 


17.6 


18 


GaN 


741 (1.0) d 


555 (0.6) d 


18.8 


18.1 


18.5 


15.1 


Si 


520 (1.0) a 


520 (1.0) a 


12.0 


13.6 


11.3 


12 


A1P 


501 c (0.5) 


440 c (1.0) 


9.7 


9.6 


7.9 


9.4 


InN 


588 (1.0) d 


457 (0.6) d 


12.6 


10.4 


8.2 


9 


Ge 


304 (1.0) a 


304 (1.0) a 


5.2 


11.7 


9.7 


8.8 


GaAs 


292 (0.6) f 


269 (1.0) f 


4.5 


8.0 


7.4 


6.8 m 


InP 


345 (0.9) f 


304 (1.0) f 


5.6 


6.0 


5.1 


5.4 


InAs 


238 s (0.5) 


217 s (1.0) 


3.4 


5.7 


4.5 


3.8 


ZnSc 


252 (1.0) h 


203 (0.7) h 


3.5 




2.6 


l.l m 


BAs 


714 (0.4) 1 


695 (l.oy 


20.8 




19.9 


19* 


GaP 


403 (1.0)° 


367 (0.2)° 


7.7 


8.9 


8.7 


9.5* 


AlAs 


404J (0.5) 


363 j (1.0) 


7.2 


8.5 


6.8 


5* 


GaSb 


237 (0.7) a 


227 (1.0) a 


3.5 


6.0 


5.6 


4.5* 


AlSb 


340 (0.3) f 


319 (1.0) f 


5.6 


4.9 


4.9 


4* 


InSb 


189 (0.7) a 


179 (1.0) a 


2.6 


4.3 


3.6 


2.2* 


ZnS 


351 (1.0) k 


276 (0.1) k 


6.2 




2.7 


1.8* 


ZnTe 


206 (1.0) 1 


177 (0.8) 1 


2.8 




2.3 


1.0* 



a Rcfcrcncc [16]. b Refcrencc [17]. c Rcfcrcncc [18]. d Reference [26]. c Refercncc [27]. f Reference [28]. 
g Reference [29]. h Reference [30]. 'Reference [31]. ^Reference [32]. k Reference [33]. Reference [34]. 
111 Reference [35]. 
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Table II. For 7 kinds of complex materials, uj rn (I m ) of vibrational Raman modes, our Vickers hard- 
ness Hy ur , and the experimental Vickers hardness Hy Xp . The Vickers hardness Hy ao predicated 
by the microscopic model has been also given as a comparison. 

Crystals mjh) u 2 {h) <Mh) m(h) cu 5 (h) Mh) tMh) H° ur H$ ao H^ xp 
/3-Si 3 N 4 1047(1.0) a 939(1.0) a 928(0.8) a 865(0.6) a 732(1.0) a 619(0.4) a 451(0.6) a 28.8 30.3 f 30 f 
a-Si0 2 464(1. 0) b 450(0. l) b 9.8 30.6 f ll f 

Stishovite 967(0.1) c 753(1.0) c 25.3 30.4 f 17-23 s , 33 s 

a-Boron 1186(0.85) d 1123(0. 12) d 933(0. 7) d 795(1.0) d 776(0. 15) d 39.8 42 h 42 h 

B 4 C 1085(1. 0) d 1000(0. 5) d 830(0.3) d 720(0.4) d 40.5 42 h 30 s , 42-49 g 

B 6 1119(0.5) e 1034(1.0) c 902(1.0)° 44.9 44 h 38 h ,45 h 

z*-BC 2 N 1328(0.9) 1326(1.0) 1292(0.3) 1176(0.6) 1076(0.4) 930(0.4) 75.5 78 f 76 f 

Reference [19]. b Reference [14]. Reference [20]. d Reference [21]. Reference [23]. f Reference [6]. 
g Reference [5]. h Reference [25]. 
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Figure Captions 
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FIG. 1: (Color online) Hardness of diamond- like crystals as a function of diagnostic Raman frequency. The 
solid line is from Eq. (4). 
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FIG. 2: (Color online) Comparison of calculated and experimental Vickers hardness for 25 crystals. Filled 
squares, circles, and triangles denote the zinc blende, wurtzite (wz), and complex crystals, respectively. 
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